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Abstract 

We present in this paper an original approach to compute the 
electrophoretic mobility of rigid nucleo-protein complexes like nucle- 
osomes. This model allows to address theoretically the influence of 
complex position along DNA, as well as wrapped length of DNA on the 
electrophoretic mobility of the complex. The predictions of the model 
are in qualitative agreement with experimental results on mononucleo- 
somes assembled on short DNA fragments (< AOObp). Influence of ad- 
ditional experimental parameters like gel concentration, ionic strength, 
effective charges is also discussed in the framework of the model, and 
is found to be qualitatively consistent with experiments when avail- 
able. Based on the present model, we propose a simple semi-empirical 
formula describing positioning of nucleosomes as seen through elec- 
trophoresis. 
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Introduction 

Electrophoresis is one of the most powerful and widely used tech- 
nique in modern molecular biology in order to address various prop- 
erties of biological samples: molecular weight and size determination 
(DNA, proteins ...), mapping of particular protein binding sites on 
DNA (enzyme footprinting) , effective charges (protein charge ladders) 
In most applications aforementioned, there is a need of some 
a calibrated sample, the so-called "ladder", in order to quantify the 
results of any electrophoresis experiments. This allows to use these 
techniques without the precise knowledge of physical mechanisms un- 
derlying electrophoresis separation. Nevertheless, by processing this 
way, one misses additional informations that are not brought by the 
comparison with the ladder. This is the case for example for nucleo- 
protein complexes like mononucleosomes. The nucleosome is the first 
degree of organization of DNA within the chromatin of eukaryotes. It 
is made from the complexation of roughly 147 bp of DNA with an 
octamer of histone proteins (Q). It has been shown indirectly that elec- 
trophoretic mobility of mononucleosomes depends on its positioning 
along DNA (Q). Now this property is widely used to detect qualita- 
tively nucleosome repositioning due either to thermal fluctuations or 
to the action of remodeling factors (0, H) . But taken alone, these kind 
of experiments do just indicate that a change occured either on the 
conformation of nucleosome, and/or on its charge distribution, since 
electrophoresis of colloidal particles is mostly sensitive to these two 
intrinsic properties. No quantitative conclusions can be reached with 
respect to the precise position of nucleosome along DNA. Physical mod- 
eling of electrophoresis might help extracting this information from the 
experiments. 

In the early works of Pennings and collaborators about position- 
dependent electrophoretic mobility of nucleosomes (jH) , datas were in- 
terpreted using similar results obtained on short bent oligonucleotides 
(6, 7): the mobility of such molecules in a gel is strongly dependent 
on both position and angle of the bent, with same qualitative trends. 
Apart from the similarity of experimental results, the systems are not 
expected to behave exactly the same, due to the large size of the nu- 
cleosome core, roughly 10 nm in diameter, which is not present in 
the bent oligonucleotides experiments. But it is quite likely that po- 
sition selectivity arises from the same physical mechanisms, still to 
be discovered. On the theoretical side, reptation models have been 
shown to explain qualitatively some features about bent-DNA but the 
position-dependent mobility cannot be obtained in a quantitative way 
(@, HI)- Again, the large size of the nucleosome core renders the rep- 
tation mechanism difficult to apply in our case, at least for explaining 
this position-dependence of mobility. Therefore there is a real specifity 
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of large DNA-protein complexes with respect to their mobility in pure 
solution or in gels, as compared to naked bent-DNA. In the present 
study, we will focus on the nuclcosomal case. The methods devcloppcd 
in this context are currently applied to analyze position-dependent 
mobility of bent DNA in a separate study (0). 

Precise computation and description of electrophoretic mobility is 
a formidable task as soon as non-trivial geometries are considered. 
Indeed one need to solve simultaneously equations describing electro- 
static potential, flow profile and ionic species distribution (|l(J ). More- 
over, any realistic model should take into account the effect of siev- 
ing medium in which electric migration is performed. We propose in 
this work a general method to evaluate electrophoretic mobility of a 
rigid nucleo-protein complex in pure buffer or in gels through effective 
continuous electro-hydrodynamic description: mimicking the confor- 
mation of nucleosome by a set of charged beads of appropriate size 
and charge, we calculate total electrophoretic mobility similarly to the 
way friction coefficients of proteins are evaluated using beads models 
mapping the protein conformation (fill . IT2T ) . This approach has been 
already applied to study theoretically the influence of different charge 
distribution on the electrophoretic mobility of polyampholytes (|l3j). 
Moreover, it will be shown that in order to reproduce quantitatively 
the experimental results specific gel features will have also to be taken 
into account in the model. Taking all these theoretical ingredients 
together allows then to investigate the influence of different physical 
factors like gel concentration, buffer ionic strength, bead-complex con- 
formation on the electrophoretic mobility. 

In order to illustrate the benefits of such an approach, we address 
two original questions in the context of mononuclcosomes characteriza- 
tion: (i) is the position-dependent electrophoretic mobility to be seen 
in pure buffer, without any sieving medium, corresponding to the case 
of capillary electrophoresis, and (u) what is the influence of nucleosome 
geometry like the amount of DNA length wrapped around the histone 
core within the nucleosome on its electrophoretic mobility? The first 
question allows to address the role of the gel in position-dependent 
mobility. In the context of the second question, the non-canonical 
conformations of a mononucleosome are supposed to mimick different 
incomplete states of nucleosomes. As an example, it is known that 
the four different histones (H2A,H2B,H3,H4) found in canonical nu- 
cleosomes are arranged into an octamer. Two different type of partial 
association of histones leading to DNA- histones complexes can also be 
found in solution: H3-H4 tetramers, and hexamers made of one H3-H4 
tetramer and one H2A-H2B dimer. These are characterized by differ- 
ent amount of DNA wrapped around the protein core. Experimentally, 
they have different electrophoretic mobilities. Another recent example 
of interest is the case of nucleosomes made of histone variants, which 
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are found in the chromatin at some specific locations along the genome 
where a strong reg ulation of gene expression occurs (either repression 
or activation) (|14T). In the case of the variant H2A.Bbd, it is believed 
that DNA wrapped length in the nucleosome variant is of order 120 
base pairs instead of the canonical 147 base pairs (|15T). Within our 
model, it is possible to evaluate the difference in electrophoretic mo- 
bility between canonical and variant nucleosomes for the same DNA 
length. 

The paper is organized as follows. In the next section, we de- 
scribe first the general formalism to compute the electrophoretic mo- 
bility of a coarse-grained nucleosome model through continuous electro- 
hydrodynamics, and then the way to include specific gel effects in this 
model. Numerical results of such models are then presented. Finally 
applications and limitations of the present work are discussed. 



Model 

General formalism 

In the present work, we propose a coarse-grained model of a mononu- 
cleosome: its shape and total charge are approximated by a rigid set 
{i} of non-overlapping charged beads of radii <jj and net charges Zi, 
hereafter denoted as the bead-complex (cf fig. Q] a) . The net steady 
state motion of such an object under an external electric field E in 
a buffer of ionic strength / and viscosity 7/ is due to the balance be- 
tween electrostatic and hydrodynamic forces. The rigidity assumption 
amounts to neglect conformation fluctuations of the whole complex 
and especially of DNA arms. This is justified for the latter as long as 
the arms are shorter than a persistence length, i.e. the thermal rigidity 
length scale, which sets the upper limit of total DNA length (wrapped 
length and arms) to roughly 400 base pairs. The neglect of confor- 
mation fluctuations of the complex is associated to the tight wrapping 
of DNA around nucleosomes. The influence of bead-complex opening 
angle fluctuations is then addressed within our model by computing 
the mobility for various rigid conformations. 

A naive statement for a rigid object like the bead-complex would 
be that the electrostatic forces are purely driving the motion, while 
the hydrodynamic forces purely exert drags, just like in any sedimen- 
tation or centrifugation experiments. However, it is well-known that 
electrostatic forces contribute as well to the net hydrodynamic drag 
due to the presence of counterions going in reverse direction of motion 
and therefore exerting an additional drag. This is the very presence of 
co and counterions that makes the problem of calculating exactly the 
electrophoretic mobility a tedious task: full solution of the problem 
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would require to solve simultaneously Poisson equation for the electro- 
static potential, Navier-Stokes equation for the flows and ion transport 
equation for the spatial distributions of ions (fiol ). 

Under certain range of parameters, it is however possible to ob- 
tain a simple closed formula by using several assumptions. The first 
one is that the distributions of co and counterions around the bead- 
complex are equilibrium distributions. This amounts to neglect the 
so-called ion-relaxation effect which is important mainly for highly 
charged objects and high electric field (|10h . The main consequence 
of this assumption is that electrostatics is now described by classi- 
cal Poisson-Boltzmann equation. The second assumption is that the 
Debye-Huckel linear approximation for the electrostatic potential is 
valid, i.e. bead-complex is not highly charged. Finally, we assume 
that the electric field driving the motion of the bead-complex is small 
enough such that orientation and polarization effects are negligible. 
The validity of these assumptions with respect to realistic systems is 
discussed in section Applications and Limitations 

Due to the linearity of Navier-Stokes equation at low Reynolds 
number as considered in this work, each bead subjected to a force 
contributes linearly to the flow field at any given point through hy- 
drodynamic interactions. Following Long el al. (jl31 ). we identify two 
different types of force on each bead, that generate different hydro- 
dynamic contributions. The first type of force {Fi} is due to the 
rigid physical connection between neighbouring beads, and it is still 
present when electrostatic interactions are switched off. The associ- 
ated long-range hydrodynamic interaction is accurately described by 
Rotne-Prager tensor (fl6lll7l). which is the first finite- volume correction 
to Oseen tensor associated to point-like forces 

Tf = (l + ^MTg (1) 



T? = _J_ (2) 

where ry is the distance between centers of beads i and j, and I is the 
identity tensor. 

The other type of force acting on each bead is electrostatic through 
the external electric field. Within Debye-Huckel approach, this gener- 
ates a screened flow profile due to the presence of co and counterions 
in the solution (18|, IjJ 2(J 21). The tensor to be used to describe this 
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The Debye-Huckel screening length k^, 1 scales with ionic strength of 
the buffer I like k d 1 - J" 1 / 2 . 

For pure translation motion, the velocity of each bead is therefore 
given by 



= v?+£(T^.F j+ T^iE) 



0, neutral , Otti 
>-0 Q,neutral 



(5) 

(6) 
(7) 



in term of the velocity v° due to each pure force field separately. 
The friction coefficient of a single bead is £° = 6irrj<7i, and the elec- 
trophorctic mobility of each bead regardless of the presence of the 
others is simply /i° = 6 *^ a . ■ Equations [5] to [7] can be cast into a 
single equation 
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where all beads have the same velocity v$ = V = fiE for a steady 
motion, and diagonal terms in the Rotne-Prager tensors are defined as 
T pp = In the case of screened hydrodynamic interactions, the 

diagonal term is given in the Appendix A ljl8l l22l,l23l). Notice that this 
term is also the inverse of isolated bead friction coefficient. 

Using the fact that the Fj's are internal forces, the final result for 
the electrophoretic mobility for a given orientation of the bead-complex 
is 



M = 



T -T, 



(9) 



where the notation of tensors T pp and T PPei has been respectively 
simplified to and T?J. The index "||" means that tensors have 
been projected along the electric field direction. Notice also that TjT^ 

is the inverse tensor of Tj^, such that J2j ^W.ij^Ujk = ^ ik ' ^ or " 
mal result was already obtained by Long et al. in their discussion of 
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polyampholyte dynamics (|13T). However, they were mainly interested 
in the influence of different charge distributions for average confor- 
mations. In the present work, the conformation is fixed due to the 
assumed rigidity of the bead-complex, and therefore Eq. [5] can directly 
and explicitely be used to calculate the electrophoretic mobility. This 
has the further advantage to keep track of bead-complex orientation 
with respect to the electric field ip. Results of the numerical calculation 
for particular geometries of the bead-complex in the case of mononu- 
cleosome are provided and discussed in the next sections. 

Specific gel effects 

Up to this point, no effect of sieving medium, i.e. the polymeric gel 
(polyacrilamide, agarose...), has been taken into account. In this sec- 
tion, we describe three main effects due to the gel, and how they can 
be incorporated into the original model: hydrodynamic flow screening, 
constrained orientation of bead-complex in the gel and trapping. 

Hydrodynamic screening - The first effect of the gel on the 
migration of bead-complex is to screen hydrodynamic flow (|18l), as 
was originally proposed by Brinkman to describe hydrodynamics in 
porous media (|24T). This effect is straightforwardly incorporated in 
the ori gina l continuous electro-hydrodynamics. Following Long and 
Ajdari (1181 ) , tensors describing screened hydrodynamic flow either due 
to electrophoretic motion or to neutral migration in a porous medium 
are identical. In the latter case, it is given by Rotne-Prager tensor in 
Eqs.[3j the Debye screening length being replaced by the gel screening 
length k^ 1 = (^gC g ) -1 / 2 , where the gel is represented as a collection 
of beads of friction coefficients £ g and concentration c g . With this 
modification of long-ranged Rotne-Prager tensor into a short-range 
one, the electrophoretic mobility of the bead complex in a gel can be 
calculated according to Eq [9l 

Orientation - The second important effect of the gel is to constrain 
the orientation of the bead-complex during its migration, see figure 
[TJd. Indeed for an anisotropic object like a mononucleosome with finite 
length DNA arms, the migration is enhanced if the size of the complex 
in the direction perpendicular to the electric field is smaller than gel 
pore size, while it is strongly reduced in the reverse situation. Using 
the continuous electro-hydrodynamic model presented in the previous 
section, this effect can be taken into account by constraining the range 
of orientation angle while performing the orientation average. This will 
be discussed more precisely in the next section. 

Trapping - Finally, the migration of nucleosomes within a gel is 
strongly influenced by trapping events: since there is a finite bending 
angle between DNA arms leaving the core of the nucleosome, this bent 
or kink might be trapped transiently through a collision with gel fibers 
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(cf fig. Q] c), just like long naked DNAs are known to hook in U- 
shape around gel fibers for high electric fields electrophoresis (25). The 
untrapping process in the case of naked DNA is thought to occur like 
a rope on a pulley. In the present case of rigid bead-complex, the 
escape from trapped configuration is mainly achieved by rigid rotation 
around the gel fiber. The overall average motion of the complex is then 
described by alternation of two states: (i) a uniform steady motion in 
the free volume of the gel with a pure buffer velocity v, corrected by the 
hydrodynamic screening effect of the gel previously mentioned, during 
an average time Tf ree , and (ii) trapping/untrapping event of vanishing 
net velocity, during an average time T trap . As a result, the average 
velocity V in the direction of the electric field is given by 

v = —^ E (10) 

Tfree 

Similar formula has been used to describe the motion of long naked 
DNA in gel when trapping events are mainly determining the overall 
dynamics f26l I27I) . 



The average time during free motion is estimated by the mean-free 
path of the bead-complex Imfp ~ l/( 7r ^ 2c s)i with d the diameter of 
gel fiber and c g its concentration. Therefore the estimate of r/ ree is 

ntPcg ( 1 1 \ 

Tfree ~ (U) 

V 

Although collision scenario is not precisely known during trapping 
events, one might anticipate that the longest time (which is the rele- 
vant time for the mobility calculation) is associated with the rotation 
of the complex around gel fibers. This motion is driven by the electro- 
static torque T e i (ip) 7 which depends on relative orientation ip between 
complex and electric field. Introducing the rotation friction coefficient 
£fl of bead-complex, the time required to escape the trap scales as 

P<P2 dip 

Ttrap ~ 6i / r , \ (12) 



t r ei(y) 



where angles (p± and if2 are respectively the orientation of complex 
at the beginning and the end of trapping event, cf figure [1}:. It will 
be checked in the next section that the precise choice of these angles 
is not that crucial to obtain qualitative informations. Moreover, the 
estimation of free and trapping average time presented here are suffi- 
cient to address the questions of position- and geometry- dependence of 
electrophoretic mobility mentioned in the introduction. Precise formu- 
lation of trapping events is beyond the scope of this work. Simulation 
works with model gels (cubic arrangement of fibers) might help to 
unravel the details of such collision events (2i 
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Results and discussion 
Migration in pure buffer 

The application of Eq. [9] for the standard geometry and conditions as 
defined in Appendix B is shown in figure [21 For a given position of the 
bead-complex, i.e. a given length for one of the arms, the mobility os- 
cillates as function of the relative orientation of electric field and bead- 
complex. This mainly reflects the anisotropy of bead-complex friction 
coefficient (datas not shown). As it is clearly demonstrated in figure 
[21 two different positions of bead-complex (x = 1 end-position and 
x = 9 middle-position) are associated with different angular average 
mobilities and oscillation amplitudes. A striking result is that in pure 
buffer, we predict slightly larger average mobility for middle-position 
than end-position bead-complex. Checking for both the electrophoresis 
and nucleosome litterature, we did not find any experimental study of 
capillary electrophoresis of positioned mononucleosomes, and therefore 
this simple prediction has not being addressed yet. 

This result in pure buffer has to be contrasted with the well-known 
experimental results obtained many times in gel, which show precisely 
the opposite: end-position nucleosomes are faster than middle-position 
ones during native gel electrophoresis. This discrepancy comes from 
the direct influence of the gel on the migration process: the porous 
structure of the gel provides an orientation constraint such that optimal 
orientation during the migration is favored rather than uniform angular 
average. As a consequence, our model predicts under such conditions 
that end-position nucleosomes move faster than middle-position ones in 
agreement with the experiments (cf thick dashed lines in figure [2|) . For 
interpolating positions of the nucleosome, the electrophoretic mobility 
changes gradually between the extreme positions. Notice that higher 
values of electric field might lead to orientation as well, due to the 
alignment of net dipole of the nucleosomes with electric field. The 
amplitude of positioning effect on the mobility is shown in the inset of 
figure [21 For the sake of simplicity, we chosed the optimal orientation of 
nucleosomes for computing the positioning curve of mobility, therefore 
neglecting fluctuations around this optimal orientation. The precise 
non-uniform distribution of orientations might be quite sensitive to 
the gel model. A qualitative comparison of the amplitude computed 
with optimal orientation with respect to experimental amplitude under 
similar conditions shows that the predicted amplitude of positioning 
effect is much weaker. Indeed ratio between middle-position and end- 
position mobility can be as small as 0.4 for particular conditions (see 
figure [5] below for longer DNA) (fiil ). In the next subsections, we 
vary different parameters of the model in order to scan the range of 
accessible amplitudes, and check whether the continuous model is able 
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to reproduce the experimental range somehow. 
Range of model parameters 

There are at least three parameters in the model than can be tuned 
in order to match experimental conditions: the hydrodynamic screen- 
ing lengths due respectively to gel and ionic buffer, and the effective 
charge of nucleosome core. For all the results presented in this section, 
optimal orientation of nucleosomes during migration has been chosen 
in order to calculate the electrophoretic mobility. The hydrodynamic 
screening due to the gel depends mainly on gel concentration for a fixed 
composition. As it is seen in figure [3J the screening length has barely 
no influence for sizes larger than the bead-complex itself. For smaller 
screening length, the positioning effect is reduced. This comes from 
the fact that as the range of hydrodynamic interactions decreases, the 
influence of position on the electrophoretic mobility decreases as well: 
for asymptotically very short range hydrodynamic interactions (of or- 
der of DNA bead size), the hydrodynamic influence of bead-complex 
arms on the core is roughly the same whathever the respective length 
of the arms. Notice that under standard conditions defined in the Ap- 
pendix B, short range electrostatic screening is always present. The 
results of Pennings et al. (f30h can be precisely interpreted as the effect 
of hydrodynamic screening: they observed that the positioning effect 
is lost when migration is performed in glycerol, therefore when the 
hydrodynamic screening is increased (without trapping). 

Coming back to the situation of no hydrodynamic screening due to 
the gel, it is possible to observe the influence of electrostatic screening 
alone on the electrophoretic mobility. As it is demonstrated in fig- 
ure [3j the positioning effect increases strongly as the ionic strength 
of the buffer is increased or equivalently as the Debye screening length 
is decreased. Note that the conformation of the nucleosome was arti- 
ficially kept fixed under ionic strength variation, in order to address 
specifically the dynamic role of salt ions. The subtle interplay between 
electrostatics and hydrodynamics provide therefore a way of modulat- 
ing the amplitude of positioning effect in pure buffer. 

Finally, the effective net charge of the core can also be considered 
as an adjustable parameter: it is difficult to assign such a value purely 
from theoretical considerations, because this net charge depends on 
many intricated features like the state of charge of the protein octamer, 
the ionic strength of buffer and the possible counterion condensation. 
What is precisely known from experiments is that the net charge of the 
nucleosome core is negative, the DNA overcharging the basic charge 
of the protein octamer. Increasing the net charge of the nucleosome 
core reduces the positioning effect, as can be seen from figure [3jThis 
reflects the role of the core as the main driving force for migration, the 
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difference in friction contribution between middle-position and end- 
position becoming less important as the net charge is increased. 

A partial conclusion drawn from this scan of model parameters is 
that continuous electro-hydrodynamic model is able to describe the 
positioning effect on the electrophoretic mobility in a qualitative way, 
but not in a quantitative way. As it will be shown in section , trapping 
events are more likely to be responsible for the experimentally observed 
amplitude of the positioning effect. 

Bead-complex opening 

Still working within the framework of continuous electro-hydrodynamic 
model, it is possible to investigate the influence of bead-complex ge- 
ometry on the electrophoretic mobility. In particular, we vary in this 
section the opening angle of the bead-complex, i. e. the amount of DNA 
wrapped in the nucleosome core, mimicking different nucleosome con- 
formations either due to incomplete formation of the histone octamer 
or to the presence of histone variants. 

During the numerical calculation, we take into account the fact that 
a reduced DNA complexed length within the nucleosome reduces ef- 
fectively its core net charge. For the sake of simplicity, we assume that 
the net charge of nucleosome core scales linearly with the complexed 
length of DNA. The results are presented in figure [4]for different open- 
ing angles 9 as function of bead-complex orientation with respect to 
the electric field for middle-position nucleosomes. The discrete values 
of 9 were chosen such that the number of beads in the arms is always 
an even number. For each opening angle, the mobility oscillates. This 
representation highlights the different amplitudes and relative phases 
of these oscillations. Due to the gradual opening of the nucleosome 
as 9 increases from 0° (two superhelical turn of DNA) to 360° (one 
superhelical turn of DNA), the optimal orientation during migration, 
as imposed by the gel, switches between two values (90° and 0°). The 
net result for the predicted mobility at optimal orientation is shown in 
the inset of figure |4j the mobility first decreases as function of opening 
angle, and increases slightly for 9 ~ 360°. 

These results only indicate qualitative trends for the opening angle 
influence for two reasons: the first one is that the variation of nucle- 
osome core net charge might be nonlinear with respect to the DNA 
wrapped length due to the ion condensation phenomena. The sec- 
ond reason is related to histone octamer stabilization by the DNA. 
In a solution under physiological conditions made by the four differ- 
ent histones, there are mainly two populations: H3-H4 tetramer and 
H2A-H2B dimers, but almost no octamer that are not stable with- 
out DNA around it. This means that if less DNA is wrapped around 
the nucleosome core, one or two H2A-H2B dimers might be lost, and 
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therefore the net charge of the core might be changed dramatically, as 
opposed to the linear variation assumed for the calculation. One way 
of addressing these two points is to test for different decreasing relation 
Z C ore v s 9. The result of such additional calculations not shown here 
is that the mobility is still a decreasing function of opening angle for 
a large range of angles: the observations presented in this subsection 
are quite robust with respect to core net charge variation. 

Trapping 

In this section, we provide a qualitative estimate of electrophoretic 
mobility according to the trapping-untrapping scenario proposed in a 
previous section about specific gel effects. A more rigorous geometri- 
cal calculation is proposed in Appendix C of this work, by evaluating 
the trap escape time due to electrostatic torque. The results of two 
approaches are consistent with each others, and therefore are thought 
to provide a correct estimate of position-dependent mobility. 

As it was previously discussed, the leading order specific gel effect 
influencing nucleosome migration is the occurence of frequent collisions 
with gel fibers. This is described approximately by introduction of two 
characteristic times (cf Eq. I10|) : the free motion time and the trapping 
time. The former depends mainly on the gel concentration, while the 
latter is closely related to the bead-complex conformation. Within 
this context, it is straightforward to interpret the fast migration of 
end-position nucleosome as compared to middle-position nucleosomes: 
the latter adopt more likely kinked configurations and their velocity are 
therefore strongly reduced through trapping, while the former adopt 
"tadpole" -like configurations and they have smaller probability to be 
trapped. 

Similarly, different opening angle lead to different trapping time. 
The limiting values are 9 — (two superhelical turns) with almost no 
trapping due to the absence of kink in the conformation, and 9 = ir 
(1.5 superhelical turn) with high trapping time due to a 180 kink. 

In order to make a simple functional prediction for the electrophoretic 
mobility, one can expand the trapping time to the second order in nu- 
cleosome position or arm length x. Taking into account positioning 
symmetry considerations, this time is rewritten 

Ttrap x(L - x) (13) 

where L is the total length of nucleosome arms and A(9) depends 
mainly on nucleosome net charge, rotational friction and geometry 
through variable 9. In the previous equation, it is implicitely assumed 
that end-position nucleosomes are not trapped at all. As a conse- 
quence, we propose the following prediction concerning the position- 
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dependent electrophoretic mobility in gel of a mononucleosome 



where fM>(%) is the mobility in pure buffer, as calculated in previous 
sections. Since variation of mobility with position in pure buffer is 
less than 6-7 % whatever the range of realistic parameters tested in 
this study, one can use constant fi$ in a first approximation in order 
to use our prediction to interpret experimental datas. Equation [14] 
provides therefore a simple semi- empirical formula that is designed to 
rationalize experimental datas of position-dependent electrophoretic 
mobility. The application of such a formula for constant /io to the 
experimental datas of Meersseman et al. (29) on mononucleosome po- 
sitioning on twofold repeat of 5S rDNA sequence (total length=414bp) 
is shown in figure El A reasonable agreement between model and 
experimental datas is found. However, the correct interpretation of 
the fitting parameter would require additional systematic experimental 
datas. Therefore we do not pursue further the analysis of fit parame- 
ters. A clear conclusion that can nevertheless be drawn at this level of 
analysis is that trapping events are mainly determining the amplitude 
of positioning effects, and therefore they cannot be neglected in the 
theoretical interpretations of datas. 



Applications and limitations 

In this work, we presented a model that can be used in order to in- 
terpret the position-dependent electrophoretic mobility of mononuclc- 
osomes. In a first step we computed the mobility in pure buffer, and 
then we took into account specific effects associated to the migration 
in gel. The theory describing the migration in pure buffer is based on 
continuous electro-hydrodynamic description as applied to a coarse- 
grained model of nucleosome made of beads of various radii and net 
charges. One of the main simplifying assumption used in order to 
derive an explicit expression out of the compact formula Eq. [9] is the 
Debye-Huckel approximation that leads to short-ranged Rotne-Prager- 
like tensors (cf eq. |3]). At first sight, it might appear very naive to 
apply it for the electrostaticpotential around such a highly charged 



object like the nucleosome (|3l[). However it allows to hide the effect 
of complicated features like counterfoil condensation in the effective 
charges of beads in the coarse-grained model. Moreover, the exponen- 
tial decay of electrostatic potential at longer range is also expected 
from more rigorous approaches. In the results presented in previous 
sections, we used the same effective charge for each DNA-bead. An 
improvement of the model would be to take into account for inhomo- 
geneous counterion condensation on the DNA-beads forming the arms 
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and therefore non-uniform DNA-bead charges, since it is known that 
electrostatic potential along a finite size polyclectrolyte is not constant 



(|32l . |33T ). The order of magnitude of changes of positioning effect on 



electrophoretic mobility due to inhomogeneous counterion condensa- 
tion is similar to the one obtained by varying the net charge of the 
core Z core (unpublished results). 

Our numerical calculations shows that the positioning effect in pure 
buffer and in gels are opposite and different in magnitude: end-position 
nucleosomes move faster (compared to middle-position nucleosomes) in 
gels, while they are slower in pure buffer. This is mainly explained by 
the orientation of the nucleosome during gel migration as opposed to 
an orientationally-averaged migration during pure buffer electrophore- 
sis at low electric field. However the small amplitude of position- 
dependent mobility in pure buffer might be difficult to measure ex- 
perimentally. A personal interpretation of our results is that capil- 
lary electrophoresis, although not used systematically for protein-DNA 
complexes characterization, might bring new information on these sys- 
tems, because both analytical and numerical hydrodynamic models 
used to interpret the datas are becoming more precise and powerful 
(see for example |13))- 

Using the continuous electro-hydrodynamic model, it is possible to 
investigate the influence of nucleosome geometry on the value of elec- 
trophoretic mobility, as it is described in previous section. The main 
result is that the mobility of middle-position nucleosome decreases as 
the DNA length wrapped in the core is decreasing. Focusing on incom- 
plete nucleosome characterization, our model would predict that the 
fastest specie in pure buffer is the octamer, then the hexamer and fi- 
nally the tetramer. Here again the experimental results are different in 
gels: the fastest are hexamer, then octamer and finally tetramer. The 
mobility of the octamer relative to the other specie is not correctly eval- 
uated through the continuous model. This discrepancy can be partially 
interpreted as a specific effect of the gel using the trapping model de- 
fined in previous section: although the precise DNA wrapped length of 
hexamer (between 1 and 1.5 superhelical turn )and tetramer (roughly 
one superhelical turn) is not known, one can speculate that the kink 
formed by the two DNA arms is less important than for the octamer 
(90° kink or 1.75 superhelical turn). Therefore the octamer is more 
sensitive to trapping mechanism, and its mobility is further reduced. 
As a consequence the octamer will not be the fastest specie anymore in 
gel. It can be either the second fastest or the third fastest. The formcr 
situation corresponds to the experimental observation. However at the 
level of the present description, we can not discriminate between the 
two cases. However, it is clear that the role of core net charge will be 
important. 

Using the trapping model, we propose a semi- empirical formula 
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Eq. [14] in order to predict the electrophoretic mobility of positioned 
nucleosomes. The application of such a formula to experimental datas 
of Meersseman et al. seems a reasonable guess f[29h . This leads us to 
propose the following method in order to determine unknown position- 
ing within a single gel run. The idea is that for this single gel run, one 
should have in a lane two different known positions (middle and end 
position for example) on a well-known sequence in order to provide a 
ladder. Then the two-parameters formula Eq. [T3](at constant /Uo) can 
be used to get the unknown positions on another sequence in different 
lanes, provided that the DNA length on which mononucleosomes are 
constructed is the same. 

As a conclusion, the main gain of the approach presented in this 
work is that it provides a rigorous framework for the understanding 
of position-dependent electrophoretic mobility of mononucleosomes. 
Moreover, the influence of different experimental parameters can be 
qualitatively predicted. This work may serve as a guideline for more 
thorough studies of electrophoresis of rigid molecular complexes. We 
are currently developping similar models in order to investigate more 
thoroughly the dependence of electrophoretic mobility of curved DNA 
on bent angle and position (Q). 
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Appendix A: friction coefficients for screened 
hydrodynamics 

In order to calculate the friction coefficient of a single bead into a fluid, 
it is necessary to solve the flow and pressure profile around this particle. 
In the case where the hydrodynamics is screened either because of 
electrostatic screening or because of the gel, the calculation of the 
friction coefficient can still be done analytically. Although the flow 
profiles are similar in the two situations, the friction coefficient are 
different due to different pressure field (fla). The reader is referred to 



the works of Russel et al. (22) or Stigter (|23l ) for futher details on 
the derivation of the friction coefficients. The result goes as follows for 
electrostatic screening: 

15 1 

fei = 67T?7eri(l + K D (J\)/[1 + — (kdcti) 2 - — (k d cti) 3 - — (k D (Ji) 4 

with the exponential integral E\ (x) — f°° e ^ du . In the case of gel 
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screening, the friction coefficient simply reads: 

igel = 67T?7<7i ^1 + K g (X\ + ^{HgCT^ 2 j ( Ki) 

Appendix B: geometry of the bead-complex 

In this appendix, we describe the geometry of the complex that mim- 
icks electro-hydrodynamics of positioned nucleosome. This complex is 
shown in figure [T] According to structural datas available on mononu- 
cleosomes ([2D, DNA is wrapped on a superhelical path around an oc- 
tamer of histones. In the case of positioned nucleosome, DNA arms 
entering and exiting the nucleosome core are also present. Their confor- 
mations depend mainly on ionic strength of the buffer (see for example 
((35h ). In the present work, we assume for the sake of simplicity that 
DNA outside the nucleosome core is following a straight path, whose 
direction is given by the tangent path of the last bead in the complex 
core . This is justified by both the rigidity of DNA backbone (per- 
sistence length of roughly 150 base pairs) for such small non-wrapped 
lengths of nucleosomes considered (< lOObp) and by the physiological 
ionic strength that effectively screens electrostatic interactions beyond 
I to a few nanometers. 

Due to the level of description of both hydrodynamic and electro- 
static interactions in this work, the nucleosome core (histone octamer 
and 147 base pairs of DNA) is represented by a single bead with an 
effective radius i? core and effective net charge Z core . Indeed, we use 
Rotne-Prager tensors as well as Debye-Huckel approach for interac- 
tions. These expression are correct for large separations, as long as 
effective radii and charges are taken into account. Moreover in the case 
of electrostatics, subtle effects like net charge of protein and counte- 
rion condensation are taken into account by the appropriate choice of 

Zcore • 

Protruding from the central core bead, DNA arms are represented 
by two linear arrays of smaller beads. Due to the natural anisotropy of 
a base pairs of radius rj^x = lnm and height r^n = 0.34rtm, a single 
"DNA" bead embeds several base pairs. For a given number of base 
pairs Nb p , the number of beads in the two arms is given by 

„„ = <^S?Ll^ + 1 (17 ) 

It has been implicitely assumed in the previous equation that in the 
reference nucleosome 147 base pairs of DNA are exactly wrapped into 
1.75 superhelical turns. This formula allows to calculate the number of 
beads in the arms for different opening angle 9, cf figure Q] The first 
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bead of each arms is tangent to the central core bead, and is located at 
the coordinate of the last base pair of the nucleosome core (base pairs 
1 and 147). Similar results are obtained if slightly different matching 
conditions are used. 

The choice of values for the bead-complex effective parameters is 
made mainly following values used in related brownian dynamics sim- 
ulations by Beard and Schlick (|36f). In the main part of this work, 
we refer "standard geometry and conditions" to the following set of 
parameters 



N, 



bp 



2506p 



R 



= 5nm 

Zcore = 200 



R 



bead 



1.35nm — no gel screening - 



Inm 
Zbead = 8.3 
no trapping— 



(18) 



Appendix C: estimation of untrapping time 
for crossed configuration 

Although, the precise collision scenario is not known, one might an- 
ticipate that the longest time (which is the relevant time to estimate 
for the mobility calculation) is associated with the rotation of bead- 
complex around gel fibers. This rotation is driven by the electric field 
that exerts a torque on the complex. 

An estimation of this torque is simply made in the case of planar 
configuration of the complex (cf figure [5]) . The result reads 



Teh 
qE 



sin ip 



6 ( Z 

1 +tan 2 - (-^-(Rcore 



Rg + 2Rbead) + (Ll + L 2 ){Rg + Rbe 



sin — 
2 



- cos ip 



L\ + L 2 ,— 4(R core + Rbead) 2 tan 2 



cos - 
2 



T 2 - T 2 



where R g = d/2 is the gel fiber radius. The charge density of DNA is 
q = Zbead /2Rbead- The length of the two arms are L\ and 

According to Eq. [12j the trapping time is mainly determined by 
integrating the inverse of electrostatic torque between two angles <pi 
and tf2 that represent respectively the initial and final orientations of 
the bead-complex during the collision. Instead of calculating such inte- 
grals, and eventually averaging over initial and final angles, we plotted 
the inverse of the torque for various arm lengths at fixed opening angle 
8 in figure \Ejp since we are mainly interested in the way trapping time 
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changes with nucleosome position: the main result is that for any cou- 
ple of reasonable angles tpi and ip2, the trapping time increases going 
from end-position to middle-position nucleosomes since curves sit on 
top of each others without crossing. As a consequence, this simple ar- 
gument shows qualitatively that end-position nucleosome will migrate 
faster than middle-position ones in a scenario where trapping deter- 
mines the dynamics. Similarly, plotting the inverse torque for various 
opening angle at fixed arm lengths shows that trapping time increases 
with opening angle in the range 9 — [0, 7r], and therefore the mobility 
decreases with the opening angle in the same range, in qualitative ag- 
greement with the results of continuous electro-hydrodynamics model. 
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Figure Legends 

Figure [H 

Gel electrophoresis of bead-complex: (a) geometry of the bead-complex 
and definition of parameters tp the orientation between complex and 
electric field, 8 the opening angle of the complex and x the length of 
one arm of the complex; the black and grey beads represent respec- 
tively DNA and nucleosome core (DNA+histones); the electric field 
direction is shown by an arrow, (b) Illustration of orientation con- 
straint on bead-complex migration in a gel. Black squares represent 
cross-section of gel fibers, (c) Illustration of two-state motion due to 
trapping- untrapping events. 

Figure [21 

Electrophoretic mobility of bead-complex in pure buffer. Main panel: 
orientation dependence with respect to eletric field of end-position 
(squares) and middle-position (circles) bead-complex. Thick grey and 
black lines represent respectively uniform angular average mobility of 
end-position and middle-position bead-complex. Dashed thick grey 
and black lines represent respectively most favorable mobility value 
(cf orientation constraint due to the gel) for end-position and middle- 
position bead-complex. Inset: Influence of bead-complex positioning 
on relative mobility ratio for most favorable orientation, x is the num- 
ber of bead in one of the arms of bead-complex. 

Figure [3]. 

Influence of hydrodynamic screening and core particle net charge. Left 
panel: Relative mobility (middle- vs end- position) as function hydro- 
dynamic screening length k^ 1 (nm) due to the gel. Center panel: Rel- 
ative mobility as function of electrostatic screening length k^ 1 (nm). 
Right panel: Relative mobility as function bead-complex core net charge 



Figure 31. 

Influence of bead-complex opening. Main panel: Mobility in pure 
buffer for standard geometry as function of orientation with respect to 
the electric field for different opening 9 — 7r/100, tt/3.5, tt/2, ir/1.25, ir* 
1.15, 7r * 1.4, 7r * 1.9. Inset: Mobility for middle-position bead-complex 
as function of opening angle 9 for optimal orientations. 
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Figure [5], 

Experimental mobility of Meersseman et al. (l29h as function of dyad 
position x pos fitted by the prediction EqOJJ The datas were obtained 
on mononucleosomes constructed on a twofold repeat of 5S rDNA (total 
length=414bp). 

Figure [61 

Estimation of trapping time T trap . (a) Geometry considered for the ro- 
tation of bead-complex due to electrostatic torque, (b) Inverse of elec- 
trostatic torque as function of orientation. Upper panel: at fixed open- 
ing angle (8 — 7r/2), different curves correspond to different asymmetry 
of arms 40bp-60bp,30bp-70bp,20bp-80bp,10bp-90bp. Lower panel: at 
fixed asymmetry of arms 40bp-60bp, different curves correspond to 
different opening angle 8 = 7r/1.5,7r/2,7r/3, ir/10. 
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